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Polymer translocation through a nanopore in a membrane investigated theoreti- 
cally. Recent experiments on voltage-driven DNA and RNA translocations through 
a nanopore indicate that the size and geometry of the pore are important factors in 
polymer dynamics. A theoretical approach is presented which explicitly takes into 
account the effect of the nanopore length and diameter for polymer motion across the 
membrane. It is shown that the length of the pore is crucial for polymer translocation 
dynamics. The present model predicts that for realistic conditions (long nanopores 
and large external fields) there are two regimes of translocation depending on poly- 
mer size: for polymer chains larger than the pore length, the velocity of translocation 
is nearly constant, while for polymer chains smaller than the pore length the velocity 
increases with decreasing polymer size. These results agree with experimental data. 



2 



I. INTRODUCTION 

Translocation of polymers across a nanopore plays a critical role in numerous natural 
phenomena and industrial processes. Many biological phenomena, such as the motion of 
DNA and RNA molecules across nuclear pores, virus infection of cells, DNA packaging into 
viral capsids, gene swapping and protein transport through membrane channels, involve the 
motion of biopolymers across membranes.!'! In chemistry, the forced permeation of poly- 
mer molecules and electrophoresis are crucial for separations and purifications of synthetic 
as well as biological macromolecules. The motion of polymers in a confined medium is 
also technologically important in food and medicine production, in oil recovery and sep- 
aration, and in many other industrial processes. Accordingly, the mechanisms of poly- 
mer translocation have become a subject of numerous experimental!®!'!!'! and theoretical 
studies.§0@0SEl@0 

A polymer molecule moving across a nanopore faces a large entropic barrier due to the 
decrease in the number of available configurations for polymer segments. In order to over- 
come this barrier and to speed up the motion of polymers, an external field or interaction 
is needed. In recent in vitro experiments,!®!®! DNA and RNA molecules are driven 
through an a-hemolysin membrane channel with the help of an external electric field. These 
elegant experiments are based on the following simple idea. When a polymer molecule moves 
through a nanopore, the electric current in the system nearly vanishes because the polymer 
blocks the flow of free ions through the channel. Accurate recordings of current blockages 
allow the description of the dynamics of translocation of single polymer molecules. The 
principal experimental findings can be summarized as follows: (i) the ability of polymers 
to enter the nanopore depends linearly on polymer concentration and exponentially on ap- 
plied voltage;!'! (ii) there is a critical value of the external electric potential below which no 
polymer molecule can enter and move through the nanopore;!! (iii) the effective number of 
free charges on a translocating polymer is surprisingly very small in comparison with the 
number of available charges;! (iv) there are two regimes of polymer threading through the 
nanopore depending on polymer length - long polymers move across the membrane with 
nearly constant velocity, while short polymers move significantly faster;! (v) the nanopore 
length defines the boundary between short and long polymers.0 These last two experimental 
observations are the subject of the present theoretical investigations. 
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Several theoretical models^'EJOoEj have been developed in order to explain these 

experimental findings, however, with limited success. Theoretical approaches to polymer 
translocation mainly follows three directions. In one approach,000 the moving polymer 
molecule should overcome the entropic barrier, and the free energies of polymer segments 
determine the dynamics of translocation. Another approachEl focuses on the interaction 
between the polymer and the nanopore, and neglects the entropic contributions from poly- 
mer segments outside the nanopore. The last approach^ views the polymer translocation 
as the motion of a kink, which travels in the direction opposite to polymer transport. All 
these theoretical works provide a reasonable description of polymer threading through the 
nanopore for very large polymers. However, these theories are less successful in understand- 
ing the dynamics of relatively short polymers due to the fact that they view the nanopore 
as an object which can hold only one monomer (with the exception of Ref.0, as discussed 
below). The a-hemolysin membrane channel, that has been used as the nanopore in in vitro 
experimentsjlfli'ii'ii has a length of approximately 5 nm for the narrow part of the chan- 
nel and thus can hold up to 10-15 DNA or RNA monomers. The theoretical approach of 
Ambjornsson et al.0 takes into account the nanopore length and studies both the polymer 
entrance into the pore and the translocation process. However, only ideal flexible polymers 
are considered and theoretical analysis focused on the dependence of polymer translocation 
dynamics on external electric field. 

In this work the effects of the nanopore length and diameter on the threading dynamics 
of single polymer molecules are investigated. The goal is to develop the simplest theoretical 
description of translocation process which takes into account the geometry of the nanopore 
and interactions between the nanopore and polymer molecule. In the present model the 
polymer moves across the membrane as shown in Fig.l. The article is organized as follows. 
In Section II we develop a model and calculate free energies and translocation times for 
translocating polymers of different sizes. In Section III we apply our results for the de- 
scription of experimental translocations of volt age- driven DNA molecules. Our theoretical 
analysis is summarized in Section IV. 
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II. MODEL FOR POLYMER TRANSLOCATION 



Consider a polymer molecule consisting of N monomers (each of size a) , which moves from 
an upper chamber to a lower chamber through a nanopore of length I = Ma and diameter 
d = Da, as shown in Fig.l. Here we assume that as soon as the polymer enters into the 
pore, it is unlikely to come back. This assumption is justified since under experimental 
conditions the energy gained by a single monomer by moving through the nanopore is much 
larger than the thermal energy, and thus the probability to return is very small.0 It is also 
assumed that the nanopore is part of an infinite two-dimensional membrane, and there are 
no interactions between the polymer and the membrane. Let the chemical potential of the 
monomer in the upper region, in the nanopore, and in the lower region be fi\, fi2 and ^ 
respectively (see Fig.l). The potential energy change is considered to occur only across the 
nanopore. 

The simple visual analysis of polymer transport across the nanopore indicates that the 
motion of long polymers (larger than the nanopore length) is qualitatively different from 
that of short polymers, and these two cases must be considered separately. In our analysis, 
we assume that both iV and M are large, which is consistent with current experimental 
conditions JfliillH In our theoretical model the translocation process starts as soon as the 
first monomer enters the pore, and ends when the last monomer leaves the pore. Note that 
experimental translocation times are slightly different, as discussed below. 



In this case, there are three regimes of motion, as shown in Fig. 2a. In regime I, the 
leading monomer enters the nanopore from the upper region and then moves across the 
nanopore. In regime II, the leading monomer leaves the nanopore, while the end monomer 
approaches the entrance of the pore. In regime III, the end monomer goes through the 
nanopore and finally leaves it for the lower region. Assuming that the polymer segments 
inside the nanopore do not contribute to the free energy, i.e., there are no fluctuations inside 
the pore, the free energy F m of the polymer configuration in regime I with m monomers in 
the pore and (N — m) monomers in the upper region, is given by0 



A. Polymers with sizes N > M 




(1) 
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where 7^ is a parameter which describes the properties of polymers and which is equal to 
0.5, 0.69 and 1 for Gaussian, self-avoiding, and rod-like chains, respectively;!^ the subindex 
1 indicates the properties of the upper region of the system. The first term in (|1]) is an 
entropic contribution due to (N — m) free monomers in the upper region, while the second 
term represents the energy gain due to moving m monomers into the pore, and includes the 
effect of the external field and chemical potential changes. The entropic contribution term 
follows from the partition function for the polymer chain in a semi-infinite space near a hard 
wall with the end monomer anchored at the wall-BEl The chemical potential difference per 
monomer is given by A/ii = /x 2 — /xi, and we assume that the potential energy inside the 
nanopore is uniform. Note that the number of monomers in the pore m can vary between 
and M. Similarly, the translocation of the polymer in regimes II and III can be described 
by 

(l-y i )ln(JV-M-m) + (l-^)lnm + ^^, (2) 



/•/;'/' ■■ '- /.7.V 

with < m < N - M, and 



— = (l- 72 )lnm + — , (3) 

with N — M < m < N, respectively. Here 7 2 describes the properties of the polymers 
in the lower region, and the chemical potential differences are A/i 2 = /J3 — /i 2 and Afi 3 = 
/i 2 - Hi + £t 3 - A*2 = A/i! + A/i 2 . 

The transport of the polymer across the nanopore can be described by a Master equation^ 

dP m (i, t) 

— m - = u m ^xP m ^ x {i, t) + w m+1 P m+1 {i, t) - (u m + w m )P m {i, t), (4) 

where P m (i,t) is the probability of moving m monomers in regime i =1, II or III at time t. 
u m is the rate constant of adding one more monomer to the segment of m monomers already 
moved, and w m is the rate constant of removing one monomer from the segment of length 
m. These rate constants are related by detailed balance, namely, 

ln = — • 5 

w m+ i k B T 

Following Muthukumar,0 it is assumed that these rate constants are independent of m; 
however, they are different for different regimes, i.e., u m = Ui for i —I, II or III, and 
generally u\ 7^ u 2 7^ m 3 . Transforming the discrete Eq. M) into continuum Smoluchovskii 
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equation, we obtain 



dP m (i,t) _ d 
dt dm 



(6) 



The mean translocation time r can now be calculated as a sum of mean first-passage times 
in each regimejl i.e., r = T\ + r 2 + r 3 , with 
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(7) 



(8) 



(9) 



N—M 

where the corresponding expressions for free energies in different regimes are used. These 
equations can be solved numerically for any set of parameters; however, explicit analytic 
results can be obtained in some special cases. Chemical potential differences are the leading 
factors in the dynamics of translocation.QE! Then, for A/i x = Afi 2 — A,u 3 — 0, we can 
obtain exact expressions for translocation times; namely, Tj are given by 



N 2 

uiii 



(l-M/N) 2 -^ 1-(1~M/Nf 



Ti 



(10) 



r 2 = a- 



(N - Mf 
u 2 



N 2 



T 3 



« 3 (2 - y 2 ) 



i 



M/N)t* 



[1 - M/Nf 



12 



(11) 



(12) 



- 1 and 



(13) 



where a is a constant, which is equal to 1/2 and n 2 /16 for the special cases j[ = 7 2 
7i = 72 = 1/2; respectively. In the limit N ^> M, these results reduce to 

M 2 aN 2 M 2 

ZUi U2 2U3 

Thus the overall translocation time r in this limit is proportional to N 2 , in agreement with 
the corresponding results from Ref.0. For another limiting case, N ~ M, we obtain in a 
similar way 



N 2 



2«i(2-7i)' 



r 2 ~ 0, 



t 3 



N 2 
2m 3 7 2 ' 



(14) 
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In this case, the overall translocation time is also proportional to iV 2 , however, with a 
different coefficient. 

For more realistic situations, when the chemical potential differences are negative and 
the entropic terms in Eqs. (|1|J^J^) are weak in comparison with the A/ij terms, we obtain in 
regime I 



in regime II, 



and in regime III, 



C' M\A^\<1, 



T-2 



fc * T(iV - M) , (N - M)\Aii 3 \ > 1, 
(JV-M)|A/x 3 |<l, 



«2|A/i 3 | 

(N-M) 2 
2u 2 ■ 



T 3 



5& M\A» 2 \ > i, 



«3]A/Lt2 

M 2 



(15) 



(16) 



(17) 



2u3 , M|A// 2 |<1. 

For large positive chemical potential differences we can easily calculate for different regimes 



1 ( k B T V / A//! 
n - — — ^— ex P M r^ 



(18) 



1 ^VexpAiV-M)^ 



M 2 \U 2 A/i £ 



A; B T 



T"3 



(19) 



(20) 



m 3 \u3An2J c \ k B T / 
When iV 3> M, the translocation time r is governed by the dynamics in regime II, and it 
becomes proportional to the polymer length for large negative chemical potential differences, 
in agreement with experimental observations.@fl 



B. Polymers with sizes N < M 

For relatively short polymers (but recall that N 3> 1), again three regimes of translocation 
are observed, as shown in Fig. 2. The motion in regimes I and III is qualitatively similar 
to the transport of long polymers (compare Fig. 2a and Fig. 2b); however the transport in 
regime II is different, since there are no polymer segments in the upper or lower regions. 
Thus the free energy expressions in regimes I and III are the same as those given by Eqs. 



8 



(H) and (!]), respectively with, however, < m < N in both regimes. The free energy in 
regime II can be taken equal to zero because at this level of approximation we neglect the 
free energy contribution from the polymer segments fluctuating inside the nanopore. 

Calculations of translocation times can be performed in a similar fashion as was done for 
long polymers. First, for translocation time in regime II at all possible values of parameters, 
it can be easily computed, 

r 2 = {N ~ M) \ (21) 

where it is assumed that u 2 ^ u 2 because the translocation process is physically different in 
this regime for short polymers in comparison with long polymers. 

When A/ii = Afi 2 = 0, the translocation times in regimes I and III are equal to 

N 2 N 2 
Tl ~ 2^(2-71)' T3 " 2^" (22) 

Here we assume that the rate constants for long and for short polymers are the same, 
since the dynamics of translocation in these regimes are very similar for short and for long 
polymers. For more realistic situations, when Afii < and A/i 2 < 0, and the entropic terms 
in the free energy expressions are small, the calculations in regime I yield 

n ~ ^ Ml|zVil1 (23) 



N 2 
2uj 



N\A^\ < 1, 



and in regime III 



" 3 2 M2 1 ^ (24) 



N 



2us , N\Afi 2 \<l. 

For large positive chemical potential differences, translocation times are given by 

2 



n = i «p^, (25) 

Mi \uxAfixJ V. k B T J 
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C. Fluctuations inside the nanopore 



So far in our calculations of polymer translocation times we neglected the contributions 
from the fluctuations of polymer segments inside the nanopore, although these fluctuations 
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may be important. To take them into account the scaling analysis can be used to describe 
the polymer molecule inside the confined cylindrical pore.il The free energy of confined 
polymer chain is given by fceTiVb, where Nf, — l/d is the number of blobs inside the pore.l§ 
Note, however, that this approach is valid when d <^ I. The size of each blob is equal to the 
diameter of the pore, i.e., 

d = Da = ag u , (27) 

where g is the number of monomers in the blob, and the exponent v is equal to 1/2, 3/5, 
and 1 for ideal, self-avoiding, and rod-like chains, respectively. Then the maximum number 
of monomers in the pore of length I = Ma and diameter d = Da is given by 

M max = l -g = MD^- l \ (28) 

Note that for rod-like chains M max = M, while for ideal flexible chains M max = MD. 

Knowing the free energy contribution of polymer segments inside the nanopore allows 
to calculate the translocation dynamics as discussed in detail above. Consider first the dy- 
namics of the polymer molecule in regime II. The contribution from the fluctuating polymer 
segments inside the pore is always constant in this regime because the number of monomers 
inside the nanopore does not change. Then this free energy term will not affect transloca- 
tion times since they are determined by free energy differences [see Eqs. (|15D , ([16D , ([171)1 • I n 
regimes I and III the number of monomers inside the nanopore is changing, however the free 
energy difference from this confinement term is equal to ksTD~ x l v ', which for experimental 
conditions!'!! is very small in comparison with entropic and chemical potential terms, and 
can be neglected. Thus the free energy contributions from fluctuating monomers inside the 
nanopore do not change the results on translocation dynamics of the polymers (provided 
that M is replaced by M max ). 



III. COMPARISON WITH EXPERIMENTS 

The above results are well compared with the experimental findings of Ref.0, where the 
size dependence of voltage-driven single-stranded DNA molecules has been investigated. In 
the present theoretical approach, the process of translocation is assumed to start as soon 
as the leading monomer enters the nanopore and to end when the end monomer leaves the 
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nanopore. Then the translocation velocity is given by 

V = (N + M)a/r. (29) 

For realistic situations (large N and M, A/ij <C 0), the results for translocation times 
in corresponding regimes [see Eqs. (|T5|) , (|T6"D , (|TTD , fl2~l~|) and (|22|)1 can be substituted into 
equation (|2"9"|), leading to explicit expressions for the translocation velocity 

\ (A r + M)«/(^ I + ^ + ^), N < M. 

However, in the experiments of Meller et al.,0 the translocation time was measured only when 
a current passing through the nanopore dropped to a level below 65% of an open channel 
current. The authors also showed that the blockade level is proportional to the fractional 
volume of the channel occupied by the polymer. This means that in these experiments the 
translocation process started when 35% of the polymer entered into the nanopore, and ended 
when only 35% of the polymer left in the pore. Thus, in order to compare our theoretical 
predictions with experimental observations, the expressions for translocation velocity ( |30|) 
should be modified as follows 

( Ar + 0.30M)a/(2« + i^£) + e S «f), N>M, 

(n + o.3<ww + ^ + taT zzr ) ■ n < m - 

Under experimental conditions,^ the single-stranded DNA molecules behave more like rod- 
like polymers, and this fact justifies using M max = M in our description of experimental data. 
Then the expressions (|3l| ) can be used to fit the observed translocation velocities! as shown 
in Fig. 3. The present theoretical approach predicts two types of translocation depending on 
polymer size. For large polymers, larger than the nanopore length, the translocation velocity 
approaches a constant value, while for short polymers the velocity increases significantly with 
decreasing polymer length. These predictions are in excellent qualitative and quantitative 
agreement with experiments for large polymers; however, for short polymers the agreement 
is only qualitative. 

There are several reasons to explain the deviations between the presented theory and 
experimental behavior for short polymers. In our theoretical approach we used a poly- 
mer description of molecule dynamics, while for such short polynucleotides (N = 4 — 12) 
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the polymer description is probably less precise and the discrete chemical nature of the 
molecules should be taken into account. In addition, our descriptions of the nanopore ge- 
ometry and the potential changes inside the nanopore are very simplified. However, for 
short polymer molecules these factors probably influence the translocation dynamics much 
stronger than for large polymers. Despite these discrepancies, the fact that a very simple 
theoretical approach can provide a qualitative and semi-quantitative description of complex 
translocation processes is rather encouraging. It also indicates that the presented theoretical 
model correctly captures and describes the main features of translocation phenomena. 

Our theoretical approach allows us to investigate the effect of interactions between the 
nanopore and the polymer molecule. The rate constants Uj measure the degree of such 
interactions. The smaller the rate constants, the larger the attraction between the moving 
polymer chains and the nanopore. As shown in Fig.4, when only the rate constants for short 
polymers in regime II are varied, the interactions between the nanopore and the polymer 
can change the translocation dynamics significantly. The stronger the interaction, the slower 
the motion of threading polymer molecules, in agreement with intuition. 



IV. SUMMARY AND CONCLUSIONS 



A simple theoretical model of polymer translocation through the long nanopores driven 
by external electric fields is presented. The fact that we take into account the nanopore 
length and diameter allows us to describe the translocation dynamics for polymer molecules 
of different sizes. By considering in detail different regimes of moving polymers across the 
membranes, the general expressions for free energies and translocation times for polymer 
chains threading through the nanopores are derived. The presented theoretical predictions 
are compared with experimental results on volt age- driven translocations of single-stranded 
DNA molecules through the a-hemolysin protein channels.@ It is found that for experimen- 
tal conditions,^ long polymers, longer than the nanopore length, translocate with nearly 
constant velocity, while short polymers move significantly faster. The theoretical analysis 
indicates that for experiments^ with a-hemolysin protein channels the polymer fluctuations 
inside the nanopore do not affect the translocation dynamics. Presented theoretical results 
are in agreement with experimental observations. 

Although a reasonable description of polymer translocation experimentsi is provided, 
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there are many factors that have not been taking into account. We assumed in our calcu- 
lations that the external field inside the nanopore is uniform, while more realistic picture 
would incorporate a potential profile inside the nanopore,0 which can be found by taking 
into account the realistic geometry of the nanopore and all electrostatic interactions. In the 
present model, the possibility that the polymer molecule can return was neglected, which is 
a very good approximation at large external driving fields, as realized in most experiments. 
Given theoretical approach allows to consider this effect by solving the Smoluchovskii equa- 
tions (H) with different boundary conditions. We also assumed that the nanopore is very 
narrow, i.e., the effect of the nanopore diameter on free energies of the polymer segments 
outside of the pore has not been considered, although the a-hemolysin pore in principle can 
hold several monomers of DNA or RNA molecules. Probably, the simplest way to include 
this possibility into presented theoretical model, is to utilize the scaling approach.0 
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Figure Captions 

Fig. 1. A polymer molecule moves from the upper chamber to lower through a cylindrical 
nanopore of length I and diameter d. The potential energy change is considered to be only 
inside the nanopore. In polymer translocation experiments (see Refs.lilffili), / ~ 5 nm, 
d ~ 2 nm, and AV ~ 50—300 mV. 

Fig. 2. Three regimes of translocation for polymers of different sizes: a) for long polymers, 
b) for short polymers. The size of the nanopore is / = Ma. Solid filled circles indicate the 
head and the end monomer of the polymer molecule. 

Fig. 3. Translocation velocity as a function of polymer size. The length of a nanopore 
is equal to 12a, where the monomer size is given by a = 4 x 10 -4 /im. Filled squares are 
experimental observations from Ref.B obtained at external field of 120 mV. Solid lines are 
our fits with A/i! = AyU 2 = 2.5k B T, Ui = u 3 = 2.3 x 10 4 s" 1 , u 2 = 1.4 x 10 4 s" 1 , and 
u 2 = 4.5 x 10 6 s -1 . Eq. ([51]) is used to calculate theoretical curves. 

Fig. 4. Translocation velocities for different degrees of interactions between the nanopore 
and the polymers. Solid curve is the same as in Fig. 3 with u' 2 = 4.5 x 10 6 , while dashed curve 
is for the case when u' 2 = 1.0 x 10 6 , dot-dashed curve is for the case when u' 2 = 0.5 x 10 6 , 
and dotted curve is for the case when u' 2 = 0.2 x 10 6 . 




Figure 1. Slonkina and Kolomeisky 




Regime I 



Regime II 



Regime III 



Figure 2a. Slonkina and Kolomeisky 
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Regime II 



Figure 2b. Slonkina and Kolomeisky 



Regime III 



